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The thermodynamic model of ultrathin lubricant film melting confined between two 
atomically-flat solid surfaces is built using the Landau phase transition approach. Non- 
equilibrium entropy is introduced describing the part of thermal motion conditioned 
by non-equilibrium and non-homogeneous character of the thermal distribution. The 
equilibrium entropy is changed during time due to transition of non- equilibrium entropy 
to the equilibrium subsystem. For the description of a melting condition the variable of 
the excess volume (disorder parameter) is introduced which arises due to chaotization 
of a solid structure in the course of melting. This variable grows with the increase in 
internal energy at melting. The thermodynamic and shear melting is described consis- 
tently. Influence of external normal pressure is analyzed on the character of lubricant 
melting. The irregular stick-slip mode of melting, which is observed in experiments, 
is described. It is shown that with growth of shear velocity the frequency of stiction 
spikes in the irregular mode increases at first, then it decreases, and the sliding mode 
comes further characterized by the constant value of friction force. Comparison of the 
obtained results with experimental data is carried out. 
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1 Introduction 



With development of nanotechnologies the friction of two smooth solid surfaces at presence of 
thin lubricant film between them is vastly investigated in last years [I]. Experimental study 
of atomically-flat surfaces of mica, which are separated by the ultrathin layer of lubricant, 
shows that lubricant has properties of solid at certain conditions j2]. In particular, the 
interrupted motion (stick — slip) is observed inherent in a dry friction ppl. Such boundary 
mode is realized, if film of lubricating material has less than four molecular layers, and is 
explained as solidification, conditioned by the compression of walls. The subsequent jump- 
like melting takes a place, when shear stress exceeds a critical value due to the effect of 
"shear melting". 

In a general case the description of behavior of ultrathin lubricant film should be carried 
out starting from the first principles, however, such approach is complicated, because the 
different lubricants and geometry of experiment are used. Therefore the phenomenological 
models, which allow to explain the experimentally observed results, are proposed. In partic- 
ular, thermodynamic (H], mechanistic (6j[7j[8], and synergetic j9] models are used. Also, the 



studies are carried out using methods of molecular dynamics 10 . It appears that lubricant 
can provide several kinetic modes j2], between which there are transitions during the pro- 
cess of friction, stipulating the interrupted friction (2j[I]. In work (7) within the framework of 
mechanistic model which is based on equation of motion for the overhead rubbing block three 
modes of friction are revealed at the account of stochastic effects: sliding mode that is proper 
low-velocity shear, regular interrupted mode, and mode of sliding at high-velocity shear. Ex- 
istence of these modes also is confirmed by the synergetic theory with taking into account 



the deformation defect of the shear modulus 11 , and the computer experiments 1 2 3 , 12 



In work [9] within the framework of the Lorenz model for approximation of viscoelastic 
medium the approach is developed according to which the transition of ultrathin lubricant 
film from the solid-like state to the liquid-like state takes a place as a result of the thermody- 
namic and shear melting. An analytical description of these phenomena was given under the 
assumption that they are results of shear stress and strain self-organization as well as of the 
lubricant temperature. Additive noises of the indicated quantities |13|14| and correlated fluc- 



tuations of temperature 15 are taken into account. Reasons for the jump-like melting and 
hysteresis, which are observed in experiments |16l|17l|18], are considered in works |19l|20l|ll|. 
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Three stationary modes are found out: two solid-like, inherent in a dry friction, and one 
liquid-like, corresponding to sliding. It is shown that transition between two last ones takes 
a place in accordance to hysteresis of dependence between stress and strain (at the jump-like 
melting), or on the temperature of friction surfaces. 

At the same time, the traditional use of the Lorenz equation set for the pointed above 
problem leads to a number of logical contradictions already on the stage of formulation of 
problem. Consideration of strain and stress as independent quantities in the meaning, that 
for each of them the separate evolutional equation is written, contradicts to the canons of 
"classic" mechanics and thermodynamics. In addition, symmetry of types of thermodynamic 
flows is absent at such formulation which predetermines strict accordance of signs in the 
mixed terms in evolution equations. An output can be found using a multidimensional ther- 
modynamic potential from which the set of Landau-like evolution equations must follow by 
standard procedure of differentiation [2T1. Earlier this approach was applied for description 



of the processes of severe plastic deformation (SPD) (22 23 24] and fracture of quasibrittle 
solids j25). Both the processes of SPD and processes ultrathin lubricant sliding have a lot in 
common, that allows supposing about legitimacy of applications of such technique in both 
cases. However, principal differences between these processes are present, mainly, related to 
an ultrathin thickness of lubricant layer (order of atomic), which bring in the limitations and 
deviations from standard procedure, which we consider in the proper place of the article. 

The general thermodynamic model of melting of ultrathin lubricant film is built in the 
offered work. Kinetic equations are written in as Landau-Khalatnikov ones for basic param- 
eters (section 2). In section 3 influence of shear velocity is considered and it is shown that 
critical value of velocity exists, at which lubricant melts in accordance to the mechanism of 
the shear melting. Influence of temperature is investigated and it is shown that at exceeding 
its critical value lubricant can melts even at zero applied shear stresses and zero velocity 
of shear, i.e., the thermodynamic melting takes place (section 4). The study of influence of 
external normal load shows that friction force grows at the increase in pressure (section 5). 
Influence of fluctuations of strain is analyzed which arise up due to errors in the experimental 
setup (section 6). All of the found features coincide with experimental data. 
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2 Thermodynamic model 



At the construction of model within the framework of Landau thermodynamic theory of 
phase transitions [21] at first it is necessary to choose a parameter the values of which char- 
acterize the examined phases. This parameter is called an order parameter and describes a 
change of state symmetry at phase transition point. An order parameter changes jump-likely 
during the first-order phase transitions and at the second-order phase transitions it changes 
continuously. However, the state symmetry changes by jump at phase transition point in 
both cases. A phase becomes more ordered with growth of order parameter, symmetry falls 
down here. 

Melting of thin lubricant film unlike melting of volume medium can take place according 
to the scenario of second-order phase transition jH]. However, there is a certain problem 
at description of the states of thin lubricants films, because such films demonstrate more 
than one type of transition p], and their states are not true thermodynamic phases, but 
they are interpreted as kinetic modes of friction III. Therefore one talks not about solid 
and liquid, but about solid-like and liquid-like phases. The increase in volume [10] and 
diffusion coefficient [l0j[27j[26] of such lubricants shows the melting. Since experimentally 
observable is the volume, for description of the state of lubricant a parameter is introduced, 
that is additional to the order parameter (the parameter of disorder) /, which has physical 
meaning of the excess volume, arising up due to the chaotization (the amorphization) of 
structure of solid during the process of melting. Unlike the order parameter usually applied 
for description of transition processes from a solid to a liquid, this parameter grows with the 
increase in internal energy during such transition. The excess volume accepts zero value at 
zero Kelvin, when all atoms of the system are densely packed and repose. It has a different 
from zero value both in the solid and in the liquid state at the real non-zero temperatures. 
However, it has a larger value in the liquid state. We introduce two values of this parameter: 
at / > fug lubricant is liquid-like, and when / < f so i, it solidifies, and the symmetry of state 
decreases here. 

Now in accordance with general procedure, it is necessary to write down expansion of 
energy in the independent parameters. Internal energy for a model, in which simultaneously 
as contributions of large shear strain so of entropy s and nonequilibrium entropy s are, 
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looks like 25 



u = M* + (T 4 + ^A(4) 2 + /i(4) 2 + as 2 + /?s4 + 

+ t ~s-h 1 S 2 + ^(Vf) 2 + V of-^if 2 -m 1 ~sf, (1) 
where Uq, cxo, A, /i, a, (3, t , ti, c, (p , <pi, m x are constants of expansion. Thus, we have: 

to = ^ + X 3 (0) 4 + Xf(4) 2 + a, S , 

h = tl+X^e^ (2) 

¥>i = ^ + 2< + /i(4) 2 ; (3) 
mi = ml+X 3 4 + ^6(4) 2 . (4) 

Elastic stresses are taken into account with accuracy to quadratic contributions via first two 

invariants of strain tensor 4? ( £ tj) 2 = e tj e % where summing is implied in repetitive indexes. 

Thus, the first invariant represents the trace of strain tensor = e\ + e| + e\, the second 

one is determined by expression [28] 

(4) 2 = (4) 2 - 2/ a = (4 + ^ + ^) 2 - 2(4^2 + e; e | + ^3) = (^i) 2 + (4) 2 + (4) 2 - (5) 

These determinations of invariants suppose that symmetric tensor 4 °f elastic strain is 
transformed to the diagonal form. 

A new basic quantity the non-equilibrium entropy s is introduced here also describing part 
of thermomotion, which is conditioned by non-equilibrium and non-homogeneous character 
of the thermal distribution. Exactly this part of entropy evolves during the process of 
external influence, tending a some stationary value. Equilibrium entropy does not evolve in 
the ordinary understanding, but changes during time due to relaxation of non-equilibrium 
entropy and its transition into an equilibrium subsystem. 

In the first line of expression (jl]) all equilibrium variables are taken into account in the 
quadratic approximation. The second line contains terms depending on non-equilibrium 
entropy, parameter of the excess volume, and equilibrium parameters. The odd degrees 
of non-equilibrium variables are positive, even ones are negative, all signs in decodings of 
coefficients are positive. 

We write down now the corresponding evolutional equations for the non-equilibrium 
parameters of state Xj in a form 

du 

T Xi Xi = ^ (6) 
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where r Xi are the times of relaxation. It should be noted that cross members in this rela- 
tionship ensue from expansion of internal energy in power series in all parameters of problem 
([IJj-Q. Thus, the coefficients of expansion together with general relaxation time in ^ give 
different kinetic coefficients for each form of non-equilibrium current. 

At description by equations ^ the system tends not to a minimum of internal energy, 
but to its maximum, that corresponds to strongly non-equilibrium processes, occurring in 
the open systems at energy pumping in them. This property of internal energy is similar 
to property of thermodynamic potential, introduced earlier for strongly non-equilibrium 



processes 29 . In our case, energy pumping is realized due to deformation at the shear of 
friction surfaces. Thus, evolutional equation for the excess volume assumes the form 

df 

T f~0^ = cV 2 f + ip - Vif - m^, (7) 

and for the non-equilibrium entropy s we obtain 

3s 

T 8 Qj; = to-tx3-m l f, (8) 

where the terms with the sign "+" describe the increase in non-equilibrium entropy due to 
the external source of energy (the work), the terms with the sign "-" reflect its going away 
to the equilibrium subsystem. 

It is needed to write down equation of evolution for equilibrium entropy. It differs from 
the usual form ([6]), since a change of equilibrium entropy occurs due to transition of its 
non-equilibrium form to equilibrium one. The decrease in non-equilibrium entropy is taken 
into account by negative terms in evolutional equation ([8]), it means that the same terms 
must take into account the increase in equilibrium entropy. Therefore evolution equation of 
equilibrium entropy is written down in the form: 

3s 

T s —=tis + m 1 f + T s L 1 (9) 

where L is the decrease in equilibrium entropy due to heat conductivity. The heat related 
to L goes away to surrounding space of lubricant, in particular, to the friction surfaces that 
play the role of thermostat. Since a decrease in non-equilibrium entropy and its transition 
to the equilibrium subsystem is the same process, the times of relaxation in expressions ^ 
and <^ are equal. Thus, non-equilibrium entropy decreases with velocity that is equal to 
rate of its transition to the equilibrium subsystem. 
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It is easily to define the current temperature of lubricant using found value of equilibrium 
entropy. According to expression for internal energy, the equation of state is valid in the 
form: 

du 

— =T = 2as + a t s + aJ + (3e e ii . (10) 

We find analytical expression for the component L in equation ([9]). In the case of non- 
homogeneous heating the equation of heat conductivity represents the ordinary equation of 



continuity 30 



where coefficient of heat conductivity k is the constant. Supposing that a layer of lubricant 
and atomically-flat surfaces have different temperatures T and T e accordingly, for a normal 
constituent V 2 Z it is possible to use the approximation with sufficient accuracy kV 2 T « 



(K,/h 2 )(T e — T), where h is a thickness of lubricant. Taking into account this, equation (11) 
is written down in more simple form 



ds K (T e \ K 

~di = K>\Y ~ l ) + T 



1 ) + -V^T, (12) 



where quantity h 2 / k plays the role of relaxation time during which smoothing of temperatures 
occurs over the thickness of lubricant due to usual heat conductivity. Accepting in (|9| 



L = ds/dt from (12), we obtain final expression for the evolution of equilibrium entropy 



Ts^ = hs + mJ + Ts-l^-lj+Ts — VlyS + T^Vlj, (13) 



where T is fixed by expression (10). According to Q, elastic stresses are determined as 

afj = du/detf 



= aoSy + Xe^ + 2iis% + (3s5 tJ + (x 3 ( % + 2X^e^ ~s - l -Xf ] s 2 5 l3 + 

+ (g6 tJ + 2^4) / - (eSij + M 3 ) f ~ (Xrfi, + 2X^8%) sf. (14) 



Expression (14) can be presented as an effective Hooke law 

a-j = aJij + 2fi eff el j + Ae-^- (15) 
with an effective elastic parameter 

/i e// = /i + X^S + flf- ^f 2 - X 6 sf. (16) 
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Constants /i e // and A are Lame coefficients (30] . The term in (15) appears that does not 
depend on strain 

1 



a v = a + Ps + Xf 5 - ^X^~s 2 + gf - ef 2 - X 3 ~sf. 



'(0) ~2 



(17) 



The first and second invariants are determined as 



n — o r . 



A + /i 



\ i] I — 133% 



Me// 



(18) 
(19) 



where n, r are the normal and tangential components of stresses acting on lubricant from 
the side of the rubbing surfaces [j] Relationships (18) and (19) represent ordinary connection 
between the components of tensors and their invariants of linear elasticity theory (see, for 
example, (28]). Here from stresses n a contribution a v is subtracted, that is not connected 
with strain. Let us use the Debye approximation relating elastic strain £?• with plastic one 



-pi. 



Total strain in a layer is determined as 



■pi 



p.. = e e. i-p* 



(20) 



f21) 



This strain fixes the motion velocity of overhead block V according to relationship 

V l' : u -''{'[, ■ i,)- (22) 



Relaxation time of strain in (20) depends on the state of lubricant: 



r e = K(j -7i/), 



(23) 



where the constants 70,71 and coefficient K are introduced. For the solid-like state of 
lubricant K = K so i. 

In the solid-like state t £ is large and therefore efj is large in accordance to the last 
expression. For the liquid-like state r £ decreases, and efj decreases. Combining relationships 



(20) - (23) we obtain expression for elastic shear strain: 

4 + 



V 



Kino - 71/) h ' 



(24) 



1 Shear stress r is defined from expression (15 1 at i ^ j, i.e., = 



Experimental data asserts that in the liquid-like state elastic strains relax very rapidly (2], 



i.e., time of relaxation for the liquid-like state is substantially smaller. Expression (23) at 



K = K so i already reflects itself the tendency of decrease in relaxation time with melting 
(at increase in /), but such dependence is fulfilled only for the solid-like state and near- 
by a transition point [H]. Therefore it is necessary to suppose for the liquid-like lubricant 

K = K liq < K sol . 

It is known that melting of lubricant has hysteresis character in most cases jl6j[T7j[T8j[7] . 
For theoretical description of the hysteresis phenomena the row of works was undertaken, 



in particular, within the framework of Lorenz model 19 , 20 , 11 . In this approach for the 



account of the indicated phenomena it is necessary to select two characteristic values of free 
volume: at / > fu q lubricant is liquid-like, and when / < f soi it solidifies. 

Let us obtain expression for friction force that is measured in experiments [2]. Besides 
elastic stresses afj the viscous ones a^ sc arise up also in lubricant. Total stress in a layer is 
the sum of these two contributions 



(25) 



Total friction force is determined by standard manner: 



F = a^A, 



(26) 



where A is an area of contact. Viscous stresses in a layer are given by the expression 31 

Msc _ VeffV 



a 



(27) 



where r) e ff is the effective viscosity that is defined only experimentally, and for the boundary 
mode 



Veff ~ {£ 



(28) 



thus for most systems 7 = 2/3. Taking into account (22), (28) the expression for viscous 



stresses (27) is written down in the form: 



a 



7+1 



(29) 



Putting ( |25| ), ( |29| in ( |26| ), we obtain final expression for friction force: 



/ IT. 



(30) 
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where of- is fixed by the expression (15) at i ^ j. According to (30), at small velocities of 
shear the lubricant is solid-like, and value afj is large. Thus, friction force takes on a largest 
value, and the dry friction is realized. At increase in velocity the interrupted mode takes 
place, when lubricant melts and solidifies periodically, accordingly friction force changes 
periodically due to the change of elastic stresses. At the further increase in velocity the 
lubricant melts, and elastic shear stresses decrease substantially, here friction force goes 
down sharply. With the further increase in velocity the value of F increases due to viscous 
component of friction force, because it becomes dominated. Thus, friction force has the 
smallest value at the critical velocity V = V c , at that the stick — slip mode disappears, and 
lubricant becomes liquid-like. Further the range of velocities is examined less or very little 
more than critical one, on which the elastic component of friction force is dominated that 
allows to study its evolution using time dependencies of elastic shear stresses. 



3 Influence of velocity and shear melting 



Ultrathin lubricant films behave differently from volume medium, therefore at their descrip- 
tion it is impossible to use standard formalism, as a row of principle new effects appears 
which must be taken into account. One of them is interrupted motion (stick — slip) (2j|4 
schematically shown in fig. 1. 
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Figure 1: Schematic realization of the interrupted (stick — slip) friction mode (2J. 

At the beginning lubricant is solid-like (stick), after exceeding some critical value of 
elastic shear stress afp it rapidly transforms in the liquid-like phase (slip) due to disordering. 
Overhead rubbing surface rises, because a change of lubricant volume takes place. The 
relaxation of cr?- occurs in the liquid-like state and lubricant solidifies again (stick) due to 
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the compressing by the walls under the action of loading. This process is periodic. One of 
the basic differences from the behavior of volume lubricants in this mechanism consists of 
that action of shear stress afj results not only in a shear but also to the increase in lubricant 
volume. This fact comports with results, which are obtained using the methods of molecular 



dynamics 10 , and can be reflected by modifying expression (18) as follows: 



The dimensionless tensor constant is introduced here, which fixes the dilatation power 
(the expansion of lubricant at a shear under action of a - ?). It is thus necessary also to 
take into account that the action of shear stresses leads to increase in lubricant thickness h. 
Relative increase in volume^] due to increase in lubricant thickness h it is possible to express 
thus: 

6V _ A6h _ 5h 

V ~ Ah ~ h ' {6 ] 

where A is the area of contact. Equating a contribution to the relative increase in volume 



from (31) due to shear stresses and the last expression, we get the change of lubricant 



thickness in the form 



Sh = h - 13 13 . (33) 



A + fj, eff 



In subsequent calculations thickness h in (24) it is necessary to replace by expression h + 8h. 
Now the model is complete, because along with the thermodynamic melting we take into 
account melting by shear. Further, within the framework of this work, for simplicity a 
homogeneous variant is examined, and V 2 = is supposed in equations Q, ( Jl3| . 

Result of cooperative numeral solution of equations Q, ([8|, (13), and p4| ) is shown in 
fig. 2. 

At zero velocity shear stress is equal to zero, the excess volume / decreases, lubricant here 
solidifies slowly due to the compressing by walls. When the system begins motion (Vi ^ 0), 
lubricant begins to melt under the action of growing stresses of,-, and the excess volume 
increases here. When / arrives at a value / = fu q , lubricant melts totaly, and as relaxation 



time in (24) becomes much smaller, stresses begin to relax. Here lubricant begins to solidify 



again, because it is supported by elastic stresses in the molten state. When it solidifies 



totally (/ = fsol), due to the increase in relaxation time in (24) the parameter / increases 



Physical meaning of the first invariant (311 is the relative change of volume SV/Vq, where SV is the 



change of volume, and Vq is the initial volume before deformation. 
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Figure 2: The time dependencies of excess volume / and shear stresses afj at parameters 
(T = 0.1 Pa, X = 200 Pa, /I = 10 8 Pa, a = 4.1 J~ x -m z -K 2 , (3 = 2.2 K,X® = 2.02 K,X% = 
4 • 10 6 K, X 4 ° = 50 J-^m^K 2 , g = 0.1 Pa, fx = 4 • 10 5 Pa, ip\ = 500 J-m" 3 , e = 10.8 Pa, p, = 
24 Pa,m\ = 3 K,X 3 = 1.12 ir,X 6 = 4040 AT, 77 = 1 J-m- 3 -s,r s = 0.1 J^-rrf-s-K 2 ,k = 
10~ n Wt-m- x -K- x ,a t = 0.2 J~ x -m i -K 2 , a v = 0.2 K,h = 10" 9 m,T e = 300 = l,n = 

-10 5 Pa, 7o = 0.0025 s, 71 = 0.005 s, K sol = 1, K Uq = 0.07, f sol = 0.04, f lig = 0.05, t* = <p* = 
0. Values of shear velocities are V — nm/s, Vi = 2 nm/s, V 2 = 4 nm/s, V 3 = 21 nm/s, 
V4 = 22 nm/s. 



again, while it does not attain value fu q , and a process repeats oneself again. According to 
this, the periodic interrupted (stick — slip) mode of melting/solidification sets in. It should 
be noted that at V = Vi the excess volume at once begins to decrease at exceeding the 
value fii q , and at solidification and achievement / = f so i it decreases some time yet, and 
it only increases after. It is related to that for an increase in / some minimum value of 



stresses is needed, and because velocity is small, it value, according to (24), increases slowly. 



Therefore, after solidification the excess volume can decrease during some time, while the 
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proper value of stresses is not attained. At increase in velocity to the value V = V 2 frequency 
of stiction spikes increases due to that at this velocity stresses rapidly increase in the system. 
Accordingly, lubricant melts rapidly, and for the same interval of time the system has time 
for accomplishing the more amount of transitions from melting to solidification. Frequency 
of peaks decreases again with the more increase in velocity V = V 3 . It takes a place due 



to that at high velocity in equation (24) the stresses relax to the greater stationary value 



at which lubricant slower solidifies. Thus, the dependence <J^(t) has a long kinetic sections 
with erf- = const. The excess volume increases during some time in this mode at exceeding 
/ > fliq and then it begins to decrease. At further increase in shear velocity V — V4 the 
interrupted mode disappears, and the kinetic mode of friction of the liquid-like lubricant 
sets in with the value of friction force F k . It takes place because with the value of velocity 
more than critical one V > V c in lubricant stress afj arises that sufficient for providing value 
/ > fsoi at which lubricant can not solidify. Let us note that at increase in velocity the 
values of stresses grow on dependence of ■(£) corresponding to the kinetic mode of friction 
with force F k . This fact comports with offered before by a mechanistic model (6). Thus, at 
the increase in velocity at first frequency of stiction spikes increases, then decreases due to 
appearance of long kinetic sections, and at exceeding the critical value of velocity V > V c 
the mode of stick — slip disappears. The described behavior well agrees with experimental 
results shown in fig. 3. 



F, mN 




v 4 >v c 



t,s 



1000 



Figure 3: The time evolution of friction force F for the hexadecane layer at temperature 17°C 
(T < T c ) with increase in sliding velocity V. Above value of critical velocity V c « 400 nm/s 
the kinetic sliding mode is realized |2|. 
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4 Influence of temperature 



Ultrathin lubricant films melt not only due to the shear melting at the increase in velocity, 
but also with increase in temperature as usual. We investigate influence of temperature on 
the examined system. For this purpose we obtain time dependencies for the excess volume 
and stresses which are similar to the shown in fig. 2. Here the value of shear velocity V 
is supposed to be a constant, and temperature of the moved surfaces T e is increased. The 
indicated dependencies are depicted in fig. 4. 

0.06 



0.04 




0.015 



0.045 



0.075 



Figure 4: The time dependencies of excess volume / and shear stresses of- at parameters 
of fig. 2 and shear velocity V = 27 nm/s. Values of temperatures of friction surfaces are 
T e0 = 50 K, T e i = 230 K, T e2 = 400 K, T e3 = 490 K, T e4 = 550 K. 

It can be seen from the figure that at the low temperatures of friction surfaces T e = T e0 
frequency of stick — slip transitions is large, and the dependence cfjit) has not the kinetic 
sections. It implies that lubricant begins to solidify after melting at once. With the increase 
in temperature (T e = T el ) frequency of peaks becomes smaller, as well as their height. 
Frequency becomes smaller due to appearance of kinetic section, i.e., lubricant now solidifies 
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slower. Decrease in peaks height implies the decline of static friction force F s . With a more 
increase in T e = T e2 the kinetic section cr|- = const becomes expressed, i.e., lubricant exists 
some time in the molten state at constant stresses which are already able to support this 
state. However, due to dissipation the excess volume decreases and lubricant solidifies, and 
the stick — slip mode is realized. At T e = T e 3 the kinetic section becomes determining, 
because here lubricant is in the liquid-like state the most of the time. At T e = T e4 lubricant 
melts finally and the kinetic mode comes. It is important to note that at the increase in 
temperature in all range a diminishing tendency for F s is saved. 



5 Influence of load 

Experiments on a boundary friction are carried out at the different values of pressure along 
the normal direction to friction surfaces. It is shown that external pressure also influences 
on the system critically j2]. Time dependencies of stresses at the constant values of system 
temperature and shear velocity are depicted in fig. 5, but at the increase in the external load 
n. 

For fig. 5a the value of temperature is chosen such that T e < T c , and stick — slip mode 
is realized. It can be seen that with increase in load friction force is increased, because the 
value grows, corresponding to the static friction. If temperature is increased to the value 
T e > T c (fig. 5b), the stick — slip mode disappears, and the kinetic mode sets in course of 
time, where with increase in pressure friction force increases also. The described features 
qualitatively coincide with experimental data obtained in work p). However, influence of 
pressure is not so simply. In particular, the increase in load can result to squeezing out the 
lubricant and, as a consequence, to decreasing the number of molecular layers between the 
friction surfaces. It causes the change of the theory parameters and the establishment of 



other friction mode 27 . Therefore for the complete account of external pressure effect it is 



necessary to carry out an additional analysis which constitutes a separate problem. 



6 Influence of noise 

The deterministic case is considered above, however, in some situations fluctuations influence 
critically on the system [14]. We consider a case, when fluctuations appear due to inaccuracy 
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0.075 



Figure 5: The time dependencies of shear stresses at parameters of fig. 2, V 
18 nm/s, and values of external normal pressure n = —0.05 MPa,ni = —0.1 MPa,ri2 
-0.2 MPa,n 3 = -0.3 MPa,n 4 = -0.4 MPa: a: T e = 200 K; b: T e = 600 AT. 



of experiment in which the value of elastic strain is badly maintained in (24), and it fluctuates. 



To that end in right-hand part of (24) we add (^-correlated stochastic source (the white 
noise mathematically defined via the Wiener process [32]), which has moments 

(34) 



<£(*)> = 0, mm)=2D6(t-t') 



where D is the intensity of source. With this addition Eq. (24) has the form of Langevin 
equation: 



+ 



V 



K(io ~ lif) h + 5h 



(35) 



the solution of which within the framework of Ito presentation with account of (34) is carried 



out with the use of iteration procedure in the form 14 



p e — <=- e _i_ 



+ 



V 



#(7o -7i/) h + 5h 
16 



At + V2DAtW n . 



(36) 



For the modeling of random force W n the Box-Muller function is used 14,33 : 



W = ^-2111^ cos(2vrr 2 ), n G (0, 1], 



(37) 



where T\ and r 2 are the pseudorandom numbers. The time dependencies of stresses at the 



parameters of fig. 2 at the use of (36) are shown in fig. 6. 

2 




0.015 



0.03 



0.045 



0.075 



Figure 6: The time dependencies of shear stresses crfj at parameters of fig. 2 with consider- 
ation of noise: a: D = 5 • 10 -5 s~ x : b: D = 2 • 10 -3 



In fig. 6a intensity of noise is small and deterministic modes are realized (fig. 2b). With the 
increase in intensity of fluctuations (fig. 6b) the stochastic regime sets in at which lubricant 
spontaneously solidifies and melts. Possibility of existence of such modes was shown by the 
methods of molecular dynamics [To]. 
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7 Conclusions 



The offered theory allows to describe effects observed at the ultrathin lubricant film melting. 
Both the ordinary case of the thermodynamic melting due to the increase in temperature and 
the shear melting are considered due to disordering under the action of the applied external 
stresses. It is shown that these two processes are connected with each other and they can 
not be examined separately. For example, at the high temperature of friction surfaces shear 
melting is realized at the smaller value of stress, and at the yet more increase in temperature 
lubricant melts even at zero stresses (the thermodynamic melting is realized). 

The stick — slip mode of friction, which is observed in experiments, is considered in model 
naturally, and its reason is the rapid relaxation of stress at achievement by lubricant the 
liquid-like state. At the temperature of friction surfaces, that is not sufficient for providing 
the melting in a state of rest, at such relaxation lubricant again solidifies and it is in the 
solid-like state during time that is necessary for the arising the stress at which melting 
occurs. Influence of pressure is also considered on the process of melting, and it is shown 
that with its increase the friction force is increased. Thus, temperature effect, shear melting, 
and external pressure are taken into consideration. These are basic factors to which pay 
attention at carrying out the experiments. The obtained dependencies qualitatively coincide 
with the experimental ones. 

It is possible to obtain a quantitative coincidence, however, as in the model there are a 
lot of parameters, it is enough labor-intensive process, and for a finding the theory param- 
eters the additional experiments must be executed. Besides, in each concrete experimental 
situation the values of parameters are different since they depend on properties of lubricant, 
friction surfaces, etc. Thus, the model is quantitative, and its modifications can be used for 
the description of concrete experiments. 

The work is partially executed with the financial support of the Fundamental Researches 
State Fund of Ukraine (Grants $25/668-2007, $25/97-2008). 
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